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ABSTRACT 

The Schonberg-Chandrasekhar (SC) limit is a well-established result in the understanding 
of stellar evolution. It provides an estimate of the point at which an evolved isothermal core 
embedded in an extended envelope begins to contract. We investigate contours of constant 
fractional mass in terms of homology invariant variables U and V and find that the SC limit 
exists because the isothermal core solution does not intersect all the contours for an envelope 
with polytropic index 3. We find that this analysis also applies to similar limits in the literature 
including the inner mass limit for polytropic models of quasi-stars. Consequently, any core 
solution that does not intersect all the fractional mass contours exhibits an associated limit and 
we identify several relevant cases where this is so. We show that a composite polytrope is at a 
fractional core mass limit when its core solution touches but does not cross the contour of the 
corresponding fractional core mass. We apply this test to realistic models of helium stars and 
find that stars typically expand when their cores are near a mass limit. Furthermore, it appears 
that stars that evolve into giants have always first exceeded an SC-like limit. 
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1 INTRODUCTION 

Once the core of a main-sequence star has exhausted its 
supply of hydrogen, it ceases to produce nuclear energy 
and, in the limit of thermal equi librium, becomes isothermal. 
ISchonberg & Cha ndrasekharl d 1942b showed that, if the envelope 
is polytropic with index n — 3, then there is a maximum fractional 
mass that the core can achieve. If the core is less massive, it can 
remain isothermal while nuclear reactions continue in a surround- 
ing shell. If this mass is exceeded then the core contracts until it 
is supported by electron degeneracy pressure or helium begins to 
burn at the centre. The idealised result is sufficiently accurate that 
it has become a well-established element of the theory of the post- 
main sequence evolution of stars. It is referred to simply as the 
Schonberg-Chandrasekhar (SC) limit. 

Simil ar limi t s hav e been computed for other polytropic 
solutions. Beechi dl988l) calculated the corresponding limit for 
an isothermal core surrounded by an envelope with n = 1. 
Eggleton, Faulkner & Cannon (1998) found that, when n — 1 in 
the envelope and n — 5 in the core, a fractional mass limit ex- 
ists if the density decreases discontinuously at the core-envelope 
boundary by a factor exceeding 3. They went further to propose 
conditions on the polytropic indices of the core and envelope that 
lead to fractional mass limits. We refer to all these limits, in cluding 
the original result of Schonberg-Chandrasekhar d 19421) as SC- 
like limits. 

Previously, we found that the black hol e mass of a poly tropic 
quasi-star exhibits a robust fractional limit dBall et al]|201 ll) . We 
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have determined why this limit exists in terms of contours of frac- 
tional core mass of solutions when plotted in the space of homology 
invariant variables U and V. We have further found that all SC-like 
limits are explained by a similar approach. In this work, we present 
our analysis, which unifies SC-like limits and indicates that they 
exist in a wider range of circumstances than the handful of cases 
discussed in the literature. In Section[2] we provide a thorough ex- 
position of the relevant features of the U-V plane. Readers who 
are familiar with these details can proceed to Section[3] where we 
present our new interpretation of SC-like limits. In Section [4] we 
provide a description that captures all the SC-like limits of Section 
fj] We also consider how to determine whether a star has reached 
an SC-like limit and how the fractional mass contours constrain its 
evolution and we conclude in Section[5] 



2 THE U-V PLANE 

This work rests on the behaviour of solutions in the U-V plane so 
we begin with a review of its features. We derive the Lane-Emden 
equation (LEE) from hydrostatic equilibrium and mass conserva- 
tion, introduce homology invariant variables U and V and explain 
their physical meaning, present the homology invariant transforma- 
tion of the equation and study the behaviour of its solutions in the 
U-V plane. We hope that, by presenting concisely the details of 
the U -V plane in a context where it is usefully applied, we might 
remove its stigma as 'that gruesome tool'Q 



1 Faulkner (2005) explains that Martin Schwarzschild described the U—V 
plane as such in a referee's report in 1965. The same quote is presumably 
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2.1 The Lane-Emden Equation 

Consider the equations of mass conservation, 



dm 2 

— — = VKT p, 

dr 

and hydrostatic equilibrium, 



dp 
dr 



Gmp 



(1) 



(2) 



where r is the distance from the centre of the star, m is the mass 
within a concentric sphere of radius r, and p and p are the pressure 
and density, respectively. We make the usual polytropic assumption 
that the pressure and density are related by p = Kp 1+ ~ , where n 
is the polytropic index and if is a constant of proportionality. We 
define the dimensionless temperatur^EI 9 by p — p c n , where p c is 
the density at the centre of the star, the dimensionless radiu jf| £ by 
r = 7/£, where 



2 (n + l)K - 
" = ^G~ pl 



(3) 



and the dimensionless mass </> by m = Airr) ' p c cf>. 

We use a polytropic equation of state to approximate a fluid 
that is between the adiabatic and isothermal limits. Shallower tem- 
perature gradients correspond to larger effective polytropic indices 
and the isothermal case (zero temperature gradient) corresponds to 
n — oo. In this case, the equation of state must be approximated 
differently but the limit is well-defined when working in the U—V 
plane. 

Certain conditions inside a star correspond to certain values of 
n. In convective zones, the temperature gradient is approximately 
adiabatic, so an ideal gas without radiation has n — 3/2 and pure 
radiation has n = 3. Real stars are more radiation-dominated to- 
wards the centre and n varies between these limiting values in con- 
vection zones. In radiative zones, n depends on the opacity. It can 
be shown, for example, that for a polytropic model with uniform 
energy generation and a Kramer's opacity l aw, n ranges f rom 13/4 
for a pure ideal gas to 7 for pure radiation jHoredtll2004 Nuclear 
burning shells and ionisation regions have shallow temperature gra- 
dients and therefore large values of n. Thus, the effective polytropic 
index can vary widely within a star. 

Introducing the dimensionless mass, temperature and radius 
into equations |QJ and (O allows us to write 



and 

d0_ 

dZ~~ 



1 



(4) 



(5) 



the citation by Eggleton et al. 1 1998) of '(Schwarzschild 1965, private com- 
munication)'. 

2 This is by the analogy to an ideal gas, for which T oc p/p. 

3 The scale factor is usually denoted a. We have used 77 to avoid confusion 
with th e density jump at the core-envelope boundary, which Eggleton et al. 
1 1998) denoted a. 



By differentiating equation l|5} and substituting for d<f>/d£ from 
equation we arrive at the LEE, 



1 d 



,d0 



e </; r 1 



(6) 



in its usual form as a single second-order ordinary differential equa- 
tion. Here, we prefer to express it as two first-order equations {4]and 
O because this preserves the physical meaning of the equations and 
easily permits arbitrary boundary conditions for the inner mass and 
radius. 

Solutions of the LEE which are regular at the centre have £ c = 
(f> c — 0. The subset of solutions that extend from the centre to 
infinite radius or the first zero of are polytropes of index n. We 
refer to solutions that are regular at the centre but truncated at some 
finite radius as polytropic cores. Conversely, solutions that extend 
from a finite radius to infinity or the first zero of 9 are polytropic 
envelopes. In addition, we refer to models that match polytropic 
cores to polytropic envelopes as composite polytropes. For n < 
5 polytropes are finite in both mass and radius while for n > 5 
they are infinite in mass and radius. The case n — 5 represents the 
threshold between the two: it has a finite mass but infinite radius. 



2.2 Homology Invariant Variables 

It is known (see IChandrasekhanll939l) that, if 0(£) is a solution 
of the LEE, then #'(£') = C^9{C£), where C is an arbitrary 
constant, is also a solution. The two solutions are homologous; the 
similarity between them is called homology. By choosing variables 
that are invariant under this transformation, we can formulate the 
LEE as a single first-order equation that captures all essential be- 
haviour. We use the variables 



U = 



and 



d log m 3p 
d log r p 



V = 



dlogp 
dlog r 



Gra p 
r p' 



(7) 



(8) 



where p = 3m/47rr 3 is the mean density of the material inside 
r. Although we have defined U and V to reduce the order of the 
LEE, the corresponding physical definitions make them meaningful 
for discussions of any stellar model. We make a brief excursion to 
explain these definitions. 

The physical variables are all positive so only the first quad- 
rant (U, V > 0) of the U—V plane is of interest. The variable U is 
three times the ratio of the local density to the mean density inside 
that radius. As r — > 0, we also have p — > p and thus U — > 3. 
We expect that the density of a stellar model decreases with ra- 
dius, so p/p < 1 and hence U < 3. The variable V is related to 
the ratio of specific gravitational binding energy to specific inter- 
nal energy. At the centre, p and p are finite and m ~ ^pr 3 , so 
V — > 0. Thus, in all physical solutions that extend to r = 0, the 
centre corresponds to (U, V) = (3, 0). If an interior solution has V 
everywhere smaller than the appropriate polytrope then it behaves 
as if it has a finite point mass at its centre.|Huntlev & S aslawl ( fl975l) 
referred to similar models, integrated outwards from a finite radius, 
as loaded polytropes. If a solution has V everywhere greater than 
the polytrope then it reaches zero mass before zero radius. Such 
solutions would have a massless core with finite radius, which is 
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Let us differentiate log U and log V as they are defined for poly- 
tropes. This gives 



Figure 1. Some general features of the U—V plane. The solid lines are, 
from top to bottom, polytropes of index 1, 3, 5 and oo. The dashed line 
is a STARS model of a 1 Mq star of solar metallicity when its radius is 
1.012 Rq and its luminosity 0.974 Lq. The dotted line shows the locus of 
the critical points G s - The locus begins in the plane at (0, 4) when n = 3 
and tends to (1,2) as n increases to oo. 



unphysical. At the surface, p — ► so U — > too. On the other 
hand, Gm/r takes a finite value but p/p oc T -> so 1/ -> oo. 
All realistic models, be they polytropes, composites of a polytropic 
core and envelope, or output from a detailed calculation, must ad- 
here to these central and surface conditions in the U—V plane. They 
therefore extend from (3, 0) towards (0, oo). Fig.Q]shows this be- 
haviour for polytropes of indices 1, 3, 5 and oo. The n = 1 and 
n = 3 models extend properly to the surface. The n — 5 and oo 
polytropes do not and therefore cannot represent a real star. In ad- 
dition, we have plotted a 1 Mq model produced by the Cambridge 
STARS code to show that it also satisfies the boundary conditions 
described above. Note that we have not calibrated this model to fit 
the Sun precisely. 

For a composite polytrope, the pressure, mass and radius are 
continuous at the join . If the density is decreased by a factor a (c.f. 
Eggle ton et alj[l9 98), then U and V decrease by the same factor. 
In other words, if p -> cT 1 p then (U, V) -> q _1 ([/, V). The 
corresponding point on the U—V plane is contracted towards the 
origin by the factor a. Such a jump occurs, for example, if there is 
a discontinuity in the mean molecular weight p between the core 
and the envelope. In this case, a = p c /p e . 

Let us now return to the polytropic solutions for which we 
defined U and V in the first place. From the definitions above, 



U 



and 



d log 4> 



(9) 



[3-n(n+ 1)"V- U] 



(11) 



i_dU _ 1 
and 

^ = \l-l + U+(n+l)- 1 V]. (12) 
The ratio of these two equations yields the first-order equation 



dV_ _ _ V ( U + (n+ 1)-V- 1 
dU ~ U \U + n{n + - 3 



(13) 



in which the dependence on £ has been eliminated. We refer to 
equation d 1 3 1 > as the homologous Lane-Emden equation (HLEE). 
The SC-like limits we wish to reproduce are shared by polytropic 
models so we now explore the behaviour of these solutions in the 
plane defined by U and V. 



2.3 Topology of the Homologous Lane-Emden Equation 

The behaviour of solutions of the HLEE is described in terms of its 
critical points , where dU /dlog £ and dV/d\og £ both tend to zero. 
lHoredtl(fl987h conducted a thorough survey of the behaviour of the 
HLEE, including the full range of n from — oo to oo in linear, cylin- 
drical and spherical geometriesQ Below, we use his convention for 
naming the critical points but consider only spherical cases with 
n ^ 1. Though realistic polytropes take n in the range 3/2 to in- 
finity, we extend it to accommodate SC-like limits discussed in the 
literature for polytropic envelopes with n = 1. 

From the numerator of equation d!3t . we see that dV/dU = 
when V — or U + V/(n + 1) = 1. The former indicates that 
solutions that approach the U-axis proceed along it until they reach 
infinity or a critical point. The latter defines a straight line in the 
U -V plane alon g which solutions are locally horizontal. Following 
lFaulknedd200^) we refer to this line as the line of horizontals. Sim- 
ilarly, from the denominator, we find dU/dV = when U = or 
U + nV/(n + 1) = 3. Again, the first locus implies that solutions 
near the V-axis have trajectories that are nearly parallel to it, while 
the second gives another straight line, this time along which solu- 
tions are vertical, hereinafter referred to as the line of verticals. The 
critical points of the HLEE are located at the intersections of these 
curves. Below, we consider the stability of the critical points as n 
varies. The analysis on which the discussion is based is provided in 
the Appendix. 

The origin is the first critical point. It is a saddle with the solu- 
tions on the V-axis approaching and those on the U -axis escaping. 
Solutions near the origin move down and to the right on the U -V 
plane. There is a further critical point on each of the axes. On the 
U -axis, Us = (3, 0) is a saddle for all values of n. It is stable along 
the U -axis and unstable across it. This point coincides with the reg- 
ular centre of realistic stellar models that we discussed previously. 
Along the V-axis, V s = (0, n+1) is also a critical point. For n < 3 
it is a source and for n > 3 a saddle. The intersection of the lines 
of horizontals and verticals is the final critical point G s ■ For each 
n, G s = ( ^-rf > 2 ) ■ The character of these points varies with n. 
Their locus is shown in Fig.Q] 



d log 9 <h 4 Readers should note that the definition of V used bv lHoredj ll987l) is 

(7l + 1) dIoii = (n+1) ^' (10) smaller by a factor n + 1. 
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Figure 2. Topology of the HLEE for n = 3. The solid line is the polytrope. 
The upper and lower dotted lines are the lines of verticals and horizontals, 
respectively. The dashed lines are a selection of solutions. Below the poly- 
trope, the solutions have a non-zero mass at the centre, which is represented 
by the critical point V B = (0, 4). Conversely, solutions above the polytrope 
have zero mass at non-zero inner radius. 



The behaviour of G s and V s distinguishes the topology of so- 
lutions into three regimes. For n < 3, V s is a pure source: it is 
unstable across and along the V-axis. The point G s has U < and 
therefore does not feature in the first quadrant of the U -V plane but 
approaches the V-axis from the left as n — > 3. When n — 3, V s 
and G 3 co-incide. The point is marginally stable across the axis. For 
n > 3, V s and G a separate. Vs is now a saddle and G a a source, 
gradually moving towards its position at (1/2,3) when n — 5. 
Fig.|2]illustrates some features of the HLEE when n = 3. The lines 
of verticals and horizontals meet at G s , which has just appeared on 
the U-V plane at (0,4). 

When n — 5, which separates the cases of finite and infi- 
nite polytropes, the U-V plane takes on a particular structure, il- 
lustrated in Fig.[3] The n — 5 polytrope is a straight line from 
Us = (3, 0) to V s = (0, 6). The point G s is a centre, with solu- 
tions forming closed loops around it. The polytrope separates so- 
lutions that circulate around G s from those that go from (oo,0) 
to (0, oo ) entirely above the polytrope. These solutions have zero 
mass at non-zero inner radius but, unlike the polytrope, have a finite 
outer radius. 

As n increases further G s becomes a spiral sink. Polytropes 
start at U s and now spiral into G s (see Fig.QJ. There is an unstable 
solution that proceeds from (oo, 0) to V a above which solutions 
extend to (0, oo). As n — > oo we also find V s — > (0, oo) and, in 
the limiting case of the isothermal sphere, all solutions ultimately 
spiral into G s because they cannot lie above the unstable solution. 



3 FRACTIONAL CORE MASS CONTOURS 

Let us consider the problem of fitting a polytropic envelope to a 
core of arbitrary mass and radius. For a given n < 5, we can regard 
a given point (Uo, Vo) in the U-V plane as the interior boundary of 



Figure 3. Topology of the HLEE for n = 5. The solid line is the polytrope. 
The dotted line from V s = (6, 0) to (1, 0) is the line of horizontals. The 
other dotted line, from (0, 18/5) to (3, 0), is the line of verticals. The point 
Gs = (1/2, 3) is now distinct from V s and is a centre, as can be seen 
from the solutions, represented by the dashed lines. Below the polytrope, 
the solutions orbit G s ■ Above, they extend from a finite inner radius, where 
4> = 0, to a finite outer mass and radius. 



a corresponding polytropic envelope by integrating the LEE from 
that point to the surface. More precisely, we can take interior con- 
ditions 



(14) 



(15) 



(16) 



Co = yJ{n+l)-W Vo, 
and 

<t>o = y/(n + 1)-*U V$ 

and integrate the LEE up to the first zero of 6 where we set £ = £i , 
This point marks the surface of a polytropic envelope, at which the 
dimensionless mass co-ordinate <j>\ is the total mass of the solution, 
including the initial value <f>o- The ratio q = cj>o/cf>i is then the 
fractional mass of a core that occupies a dimensionless radius £o- 
By associating each point in the U-V plane with the value of q for 
a polytropic envelope that starts there, we define a surface q(U, V). 
We use the contours of this surface to characterise the SC limit. 

Figs|4] and [5] show contours of q(U, V) for polytropic en- 
velopes with n — 3 and n = 1 respectively, along with a selection 
of interior solutions that lead to SC-like limits. For n < 3 the con- 
tours are dominated by the critical point V s and for n > 3 by G s . 
Away from V s or G s all the contours at first curve away from the 
U -axis and then tend towards straight lines. 



3.1 The Schonberg-Chandrasekhar Limit 



iKippenhahn & Weigertl<ll990h discuss the SC limit in terms of frac- 
tional mass contours, jcannonl fl992T) also explicitly described the 
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Figure 4. The dashed lines are contours of a core's fractional mass q = 
4>0 /<t>l beneath an envelope with n = 3. They increase in steps of 0. 1 from 
0.1, at the bottom, to 0.9, at the top. The larger solid spiral is the isothermal 
core with a = 1. The smaller spiral represents an isothermal core when 
a = 2. The upper and lower straight lines represent the inner boundaries 
for quasi-stars with b = 3 and b = 1 respectively (see Section l33l . 



Figure 5. The dashed lines are contours of fractional mass q = ipo/<f>i for 
n = 1. The solid spiral is again the isothermal core solution. The top-most 
diagonal line is the polytrope of index 5 with a = 1. The other diagonal 
lines are, from top to bottom, core-envelope bound ary conditions for the 
envelope when a = 2,3,4 for n = 5 as shown bv lEggletonetalJh998l 
Fig. 3). 



SC limit in terms of fractional mass contours, although he em- 
ployed a different set of homology-invariant variables. Fig.[4]shows 
the isothermal solution along with the fractional mass contours for 
n = 3 envelopes. The SC limit exists because the isothermal so- 
lution only intersects fractional mass contours up to a maximum 
g m ax = 0.359 when a = 1. In other words, along the isothermal 
solution, the function q(U, V) achieves a maximum of 0.359. 

A maximum V for the isothermal curve exists because, for a 
given mass, there is a finite maximum pressure that a core can ex- 
ert. This limit is usually derived by defining the core pressure using 
viria l arguments and maximizing it with respect to the core radius 
(e.g. lKippenhahn & Weigert 1990). Such an explanation partly de- 
scribes the SC limit but our interpretation makes clear that the ex- 
istence of the SC limit has as much to do with the behaviour of the 
envelope solutions as the isothermal core. For example, changing 
the polytropic index of the envelope changes the mass limit. 

If there is a density jump by a factor a at the core-envelope 
boundary (see Section 2.2), then U and V at the edge of the core 
must be transformed to find the base of the envelope in the U—V 
plane. That is, if p -> a" 1 p, then (U, V) -» a _1 (f/, V). The 
contraction of the isothermal core for a = 2 is included in Fig.|4] 
The inner boundary of the envelope shifts to a smaller fractional 
mass of about 0.09 so the SC limit falls too. 

The argument presented here implies that SC-like limits ex- 
ist whenever an envelope is matched to a core that only intersects 
fractional mass contours of that envelope up to some maximum. 
We now use this to explain the existence of other mass limits in the 
literature. 



3.2 Related Polytropic Limits 

Fig-H]shows the fractional mass contours for n = 1, the isothermal 
solution, and n = 5 polytropes w ith a — 1, 2, 3, 4 as used by 
lEggleton et aHdl998h . lBeech1 Jl988l) calculated an SC-like limit for 
an isothermal core embedded in a polytropic envelope with n = 1. 
Because the behaviour of fractional mass contours is similar for 
n = 1 and n — 3, the existence of t he limit i s now no surprise. Note 
that the numerical value found by iBeechl (1988) differs because 
he included the radiation pressure of the isothermal core, which 
displaces U and V at th e core-envelope bou ndary. 

The conclusions of lEggleton et al.l j 19981) are also catered for. 
The critical point V s = (0, n + 1) separat es solutions, and thus 
contours, with q « from those with q w 1. Eg gleton et alj ( fl998l) 
concluded that, for n = 1 envelopes, cores with n < 5 are never 
subject to an SC limit; those with n > 5 always are; and those with 
n = 5 constitute the margina l case for which the limit exists when 
a > 3. Eggleton et al. ( 1998) extended their arguments to changing 
the polytropic index of the envelope. Because these conclusions are 
based on the critical behaviour of the solutions, which is reflected in 
the behaviour of the contours, the same results follow here. We have 
shown how they are characterised by the contours in the same way 
as other limits and are a particular example of our broader result. 
That is, we have shown that SC-like limits exist whenever the core 
solution fail s to intersect all frac tional mass contours. The cases 
identified by Eggl eton et alj {l998) fall within this description. 

3.3 Loaded Polytropes and Quasi-stars 

Quasi-stars are objects consisting of a stellar-mass black hole em- 
bedded in a massive, hydrostatic, giant-like envelope, potentially 
formed when primord ial gas in large dark m atter haloes collapsed 
in the early Universe ( IBegelman et al.l2 006). The black hole is able 
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to grow rapidly as long as the hydros tatic structure persists. Follow- 
ing th e simple models de scribed bv lBegelman. Rossi & Ar mitage 
d2008h . lionet all l l201lh computed models with the Cambridge 
STARS code and found a maximum mass for the black hole that 
was accurately reproduced by polytropic models. We now show 
how this result is related to our analysis of the SC limit. 

The interior boundary condition for the quasi-star models can 
be written as 



r = ~- 



1 2Gm 



(17) 



i a fraction 1/6 o f the Bondi r adius, where mc, /2 = Gm/r. 
nan et alj J2008h used b = 3; lBall et alj feOl lh used 6=1. 



where 6 is a scale factor, mo is the mass interior to ro and 
c s = IV I P is the adiabatic sound speed. The boundary condition 
is then i 
iBegelman i 

Accretion on to the central black hole supports the envelope by ra- 
diating near the Eddington limit of the entire object so the envelope 
is strongly convective and the pressure is dominated by radiation. 
The envelope is approximately polytropic with index n = 3. Now, 
at the interior boundary, (? s — (1 + l/n)Kpc/ n , mo = 47rr/ 3 p c (j>o, 
and ro = 77^0 so 



<j>o 



2n 



€0. 



(18) 



Transforming to U and V gives U = 2nfo/& and V = b(n + 
l)/2n. Varying £0 traces a straight line, parallel to the t/-axis. In 
Fig-S we have plotted Vo = 2/3 and 2, which correspond to b = 1 
and 3, respectively, for n — 3. The line of Vo does not intersect 
all the contours of fractional core mass because many of them are 
positively curved. Thus, a mass limit exists, as in previous cases. 
For larger values of 6, Vo is also larger and intersects more of the 
contours. The mass limit is therefore larger. 

We have limited ourselves to the case where n — 3. The frac- 
tional mass contours in Figs|4] and [5] show similar behaviour. Con- 
vective envelopes are approximately adiabatic and have effective 
polytropic indices between 3/2 and 3, depending on the relative 
importance of gas and radiation pressures. All such envelopes pos- 
sess fractional mass contours that are similar to the two cases here 
and we conclude that a fractional mass limit for the black hole ex- 
ists in all realistic cases. Envelopes with 3 < n < 5 have more 
complicated fractional mass contours so we cannot immediately 

draw similar conclusions. 

This mass limit is not exactly the same as found bv lBall et al.1 

J201 ]]). They computed the Bondi radius using the mass of the 
black hole only, even once the mass of gas inside the Bondi ra- 
dius is comparable to (and even exceeds) the mass of the black 
hole. Begelman (2010, private communication) pointed out that the 
Bondi radius should be defined for the total mass inside ro, not just 
the black hole mass. Presuming that the gas has a density distri- 



bution p(r) 
gives the equation 



-3/2 



inside the cavity around the black hole, this 



ro = 



2G 



87T 3 

Mbh + -^-por 



(19) 



where Mbh is the mass of the black hole only, for b = 1. Making 
the same substitutions as in equation d 1 8h for polytropic index n = 
3, the equation becomes 



i>BH + \ £0 - \io = 0, 



(20) 



which only has a real positive root if 0bh < l/(18\/3). The corre- 
sponding fractional core mass limit is q — 0.0 166. 

In trying to move ro inwards, we found jBall et alj|20110 we 
could not construct models with b 3.8 in equation (T7J. We 
can see how this comes about from the behaviour of the polytropic 
limit. As b increases, Vo increases and eventually passes the critical 
point G s when the nature of the limit is reversed. When b > 2n, 
small inner masses correspond to envelopes with negligible enve- 
lope mass so it becomes impossible to embed a small black hole 
inside a massive envelope. In other words, the m ass limit become s 
a minimum inner mass limit. For the models of iBall et alj d201 ll) . 
the finite mass of the black hole corresponds to a finite value of Uo 
that displaces the envelope slightly from the V-axis. The fractional 
mass contours are closely packed near G s , so a small value of Uo 
introduces a minimum inner mass limit for b < 2n. 



4 GENERAL LIMITS 

The SC-like limits discussed above all exist because each locus of 
core-envelope boundaries only intersects fractional mass contours 
with q smaller than some q mwL . We can use this condition to iden- 
tify large classes of core solutions that lead to SC-like limits. For 
example, an SC-like limit must exist whenever the core solution 
has everywhere V < n e + 1, where n e < 5 is the polytropic in- 
dex of the envelope. Cores desc ribed by n = 5 poly tropes with 
a > 6/(n e + 1), as discussed by Eggleton et al. ( 1998), satisfy this 
condition. SC-like limits also exist whenever the curve defining the 
inner edge of the envelope touches but does not cross some frac- 
tional mass contour. This explains the original SC limit and also 
implies that any composite polytrope with n > 5 in the core and 
71 < 5 in the envelope is subject to an SC-like limit. In fact, the lim- 
its are determined by behaviour of solutions at the core-envelope 
boundary, so SC-like limits also exist when there is a layer with 
n ^> 5 at the base of the envelope. The constraint becomes even 
stronger as the density gradient at the core-envelope boundary be- 
comes steeper or the mean molecular weight jump becomes more 
pronounced. Both of these conditions become relevant immediately 
after a star leaves its core-burning sequence, which implies that SC- 
like limits apply earlier in a star's life than previously thought. 

For a given core solution and polytropic index of the envelope, 
there are potentially two solutions either side of an SC-like limit 
that correspond to the same fractional core mass. The evolutionary 
sequence of static models determines which solutions occur in re- 
ality. For example, quasi-stars are initially constructed by loading 
a star with a small mass q <C 10 -3 . These lie close to the V-axis 
so solutions must start with Uo « on the locus of inner boundary 
conditions and move towards greater Uo as the black hole grows 
and the quasi-star evolves. The evolution halts when the maximum 
mass is achieved because the black hole cannot lose mass. Further 
solutions, with smaller inner masses than the maximum, exist as 
Uo increases further. We have calculated polytropic models with 
Uo greater than the SC-limited value but we were not able to com- 
pute STARS models along the same sequence. 

It is also possible to test whether a given composite polytrope 
is at an SC-like limit. If the core solution touches but does not cross 
the contour corresponding to the given fractional core mass, then 
the model is at an SC-like limit. If this condition is satisfied, then 
extending or contracting the core can only admit a smaller frac- 
tional core mass. The test can be applied to realistic stellar models 
but identifying stars that have reached an SC-like limit is difficult. 
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Figure 6. The solid lines are Cambridge STARS models of a 1 Mq pure 
helium star. Evolution proceeds from the right-most line at large V to the 
left. The short-dashed lines are, from bottom to top, contours of fractional 
masses 0.516, 0.543 and 0.562 for n = 3. The inset is a magnification 
of the boxed region. The first model is at the end of the core He-burning 
sequence. The second model appears to have reached an SC-like limit of 
q 0.543 according to our criterion. At this point, the star begins to ex- 
pand rapidly and T c ff decreases. As it evolves across the Hertzsprung gap, 
the gradient dV/dU becomes steeper than the contours for the relevant 
fractional core mass. 

The inner edge of the envelope is not clearly denned and the effec- 
tive polytropic index varies throughout the envelope. 

Despite the approximate nature of the test, we have applied 
it to models produced by the Cambridge STARS code. Figs[6]and 
[7] show models of pure helium stars of IMq and 0.5 Mq. The 
envelopes are still mostly radiative and their effective polytropic 
indices vary between n = 2.5 and 3.5, so we have used fractional 
mass contours for envelopes with n — 3 everywhere. The mass 
co-ordinate at the core-envelope boundary was determined approx- 
imately by eye at the point where V was at a local minimum. Where 
a minimum did not exist a similar nearby point in the models was 
taken around where dV/dU was a maximum. Based on these, both 
stars appear to reach SC-like limits shortly after moving off their 
core-burning sequences although their subsequent evolution is dif- 
ferent (see Section l4~2l . The regions which lie parallel to the con- 
tours have n> 5. Although the n = 3 contours are an approxi- 
mation, it appears that both stars reach SC-like limits shortly after 
leaving their core-burning sequence. 

4.1 Beyond the limit 

What happens when an isothermal core exceeds an SC-like limit? 
In short, its effective polytropic index must change. This can hap- 
pen in two ways. Under suitable conditions, the inner part of the 
core becomes degenerate. Degenerate matter is described by a poly- 
tropic equation of state with n — 3orn = 3/2 in relativistic and 
non-relativistic cases, respectively. The inner core can tend to such 
an equation with an isothermal layer further from the centre. An 
SC-like limit still exists but the isothermal layer is displaced up- 
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Figure 7. The solid lines are Cambridge STARS models of a 0.5 Mq pure 
helium star. Evolution proceeds from the right-most line at large V to the 
left. The short-dashed lines are, from bottom to top, contours of fractional 
masses 0.640, 0.740 and 0.825 for n = 3. The inset is a magnification 
of the boxed region. The first model has already left the core He-burning 
sequence. The second model appears to be at an SC-like limit of g max ~ 
0.74 but, although the star is expanding, it does not evolve into a giant. The 
third model is at the maximum luminosity achieved by the model and is still 
at or near an SC-like limit. As the star continues to evolve, its profile in the 
U—V plane tends to a polytrope with n es 1.5. The star does not become a 
giant and evolves directly on to the white-dwarf cooling sequence. 

wards in V, so a larger core mass is possible. Alternatively, the 
core departs from thermal equilibrium and contracts. For an ideal 
gas, dp/dp = 1 + 1/n = 1 + dT/dp, so the temperature gradient 
decreases the effective polytropic n. As long as n > 5, an SC-like 
limit persists but, as in the previous case, it corresponds to a larger 
fractional core mass. 

The structure of the envelope offers some respite from the con- 
straints imposed by the core. For radiative envelopes, where n is 
not much greater than 3, the fractional mass contours are much like 
those shown in Fig.|4] If the envelope becomes convective, then the 
effective polytropic index varies between 3/2 and 3. An SC-like 
limit still exists but the behaviour of the fractional mass contours is 
less extreme near V 3 for smaller values of n (compare Figs|4]and|5](- 
Away from G s , contours run along smaller values of V for smaller 
n. Equivalently, q(U, V) is larger at a given point (U, V) for small 
values of n. For example, for n = 3, q(2, 4) = 0.528, whereas for 
n = 3/2, q(2, 4) = 0.562 and for n = 1, q(2, 4) = 0.576. Thus, a 
smaller polytropic index in the envelope permits a larger fractional 
core mass. 

4.2 Evolution into giants 

Although the evolution of main-sequence stars into giants is re- 
produced by detailed calculations of stellar structure, the cause 
of a star's substan tial expansion after leavin g the main sequence 
remains unknown. Eggl eton & Cannon! jl99lf) showed that if the 
effective polytropic index of a star is everywhere less than some 
fimax < 5, then it is less centrally condensed (i.e. p c /p is smaller) 
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than the polytrope of index n m ,-n . lEggletonl J2000l) further conjec- 
tured that the evolution of dwarfs into giants during shell burning 
therefore requires that a significant part deep in the envelope has an 
effective polytropic index n>5. This condition is similar to the 
condition under which a star is subject to an SC-like limit and the 
phenomena may be related. 

Fig.[6] shows the evolution of a IMq helium star. After core 
He-burning is complete, the star briefly contracts and nuclear fu- 
sion continues in the portion of the core that was not convective 
during its burning phase. The star expands slightly before appear- 
ing to reach an SC-like limit. At this point, the expansion of the 
envelope accelerates and the effective temperature decreases, indi- 
cating that it has begun to evolve into a giant. The star continues to 
move across the Hertzsprung gap and the calculation is terminated 
when C-burning begins at the centre of the core. We also evalu- 
ated the evolution of solar metallicity stars with masses 1, 3, 5, 
7 and 9 Mq and found qualitatively similar behaviour. All appear 
to reach an SC-like limit at the end of the main sequence, before 
starting to cross the Hertzprung gap and evolving into giants. 

Fig.[7]shows the evolution of a 0.5 Mq helium star. The star 
appears to reach an SC-like limit in its shell-burning phase when 
the fractional core mass is about 0.74. The star expands briefly but 
the U-V profile tends back towards a polytrope thereafter. If the 
mean molecular weight gradient were steeper then an SC-like limit 
would apply earlier because the gradients of the contours are shal- 
lower at smaller V. However, the density jump between helium and 
metals is modest and allows a large core to develop before a limit 
applies. The star's surface temperature does not decrease until it 
reaches the white dwarf cooling sequence. Thus, although this star 
appears to reach an SC-like limit, it does not become a giant. 

Both these stars appear to reach SC-like limits but only one 
becomes a giant. What is the crucial difference between them? In 
the 0.5 Mq star, the expansion around the limiting point reduces 
the effective polytropic index in the burning shells and the limit re- 
laxes slightly. In the 1 Mq star, the burning shell is thinner, hotter 
and deeper within the envelope. The shell responds less to the ex- 
pansion. It appears that reaching the limit leads to expansion of the 
envelope but the response of the structure may then release the star 
from its SC-like limit as for the 0.5 Mq star and not the 1 Mq . 

Another difference between the two stars is that, following the 
limiting point, the core of 1 Mq star extends beyond the SC-like 
limit. The gradient of the core profile becomes steeper than the 
fractional core-mass contour, unlike the 0.5 Mq star where it is 
approximately equal and then becomes shallower. Such a structure 
is not limited because, by reducing the radial size of the limiting re- 
gion in the U-V plane, a larger fractional core mass can be accom- 
modated but in order to reach this point it must have been limited 
before. 



We also derived a test of whether a polytropic model is at an 
SC-like limit. If the core solution touches but does not cross the 
contour of the appropriate fractional mass, then the model is SC- 
limited. Although the condition is only approximate for realistic 
models, we have applied it to helium stars and found that achieving 
a limit corresponds with expansion of their envelopes. Stars that 
clearly exceed an SC-like limit consistently evolve into giants but 
it is not known if this connection is causal. We have thus demon- 
strated that the original SC limit is a particular case of a broader 
phenomenon, that SC-like limits apply earlier in a star's evolution 
than previously thought, and that there is a connection between ex- 
ceeding these limits and evolving into a giant. 
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APPENDIX A: 

From equations J 1 1 b and J 1 2b . we find 



5 CONCLUSION 

We have shown that SC-like limits exist whenever the solution de- 
scribing a stellar core only intersects a fraction of contours of con- 
stant fractional mass for envelopes with a given polytropic index. 
O ur descr i ption expl ains the original SC lim it, SC-like limits found 
bv lBeechl Jl988l) andlEggleton et all dl99 j), and the limit for poly- 



tropic quasi-stars found by Ball et al.l d201 lb . It also shows that SC- 
like limits exist under a wide range of circumstances. This includes 
models where the core solution touches but does not cross a partic- 
ular fractional mass contour, as is the case if there is a layer with 
n 2> 5 at the base of the envelope. 



dU 
dloge 
and 

dV 



-U[U + n(n+ l)"V-3] 



V[U + (n+ 1)"V- 1]. 



(21) 



(22) 



dlog£ 

This is an autonomous system of equations: the derivatives depend 
only on the dependent variables U and V. The linear behaviour of 
such systems around the critical points can be characteris ed by the 
eigen vectors and eigenvalues of the Jacobian matrix (e.g.[Strogatz 
1 1994b , 
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Table 1. Critical points of the HLEE. A n = yfl + n(22 - 7n). 

Critical point Eigenvalues Eigenvectors 

Os (0,0) 3 -1 (1,0) (0,1) 

Us (3,0) -3 2 (1,0) (~3n, 5 + 5n) 

K (0,n + l) 1 3-n (0,1) (2-n, 1 + n) 

Gs fe,2^±i) (l-n T A„,4 + 4n) 



J= (23) 

\ at/ d log 5 dV d log 5 / 

V V -l + tf + ^V 

at the critical point in question. In particular, if the real component 
of an eigenvalue is positive or negative, solutions tend away from or 
towards that point along the corresponding eigenvector. Such points 
are sources or sinks. When the point has one positive and one neg- 
ative eigenvalue it is a saddle. If the eigenvalues have imaginary 
components then solutions orbit the point as they approach or re- 
cede. We describe these as spiral sources or sinks. If the eigenval- 
ues are purely imaginary, then solutions form closed loops around 
that point, which we call a centre. The choice of independent vari- 
able, in this case log£, is not relevant in such analysis. 

Table 1 shows the eigenvalues and eigenvectors for the criti- 
cal points in Section l231 as functions of n. The origin O s is always 
a saddle, with paths approaching along the V-axis and escaping 
along the (/-axis. Because n ^ 1 for realistic or interesting mod- 
els, Us also keeps the same saddle behaviour in our discussion, 
with points approaching along the axis and escaping along the other 
eigenvector, which always points towards the top left of the U—V 
plane. For n < 3, V s is a source. One eigenvector is always along 
the V-axis and the other across it but the latter varies from pointing 
up in the U -V plane to pointing down. When n = 3 one eigenvalue 
is zero so that it is a point of marginal stability. Points along the cor- 
responding eigenvector are also stationary in the linear regime. As 
n increases past 3, V s becomes a saddle, with points now approach- 
ing from positive U. 

Lastly, G s displays the most complicated behaviour. It first 
appears in the U-V plane when n = 3. In this case, it co- 
incides with V s . As n increases, G s is at first a source. When 
A n = yl-j- n(22 — In) = 0, the eigenvalues take on an imag- 
inary component, so G s becomes a spiral source. Increasing in n, 
the special case n — 5 is reached. The eigenvalues become purely 
imaginary at Gs, so the point is a pure centre (see Fig.[3jl. For 
n > 5, the real part of G s is negative and it becomes a spiral sink. 
As n — > oo, Vs effectively vanishes and all solutions ultimately 
reach G s . 
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